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We investigate the simultaneous estimation of n not necessarily identical unitary channels using 
multi-partite entanglement. We examine whether it is possible for the rate at which the mean square 
error decreases to be greater than that using the channels individually. For a reasonably general 
situation, in which there is no functional dependence between the channels, we show that this is 
not possible. We look at a case in which the channels are not necessarily identical but depend on 
a common variable. In this case, the mean square error decreases more rapidly using multi-partite 
entanglement. 
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Introduction. Estimation of quantum states and chan- 
nels is of fundamental importance to quantum informa- 
tion theory. Estimation of unitary channels, when n 
copies are available, has received a lot of attention; for an 
accessible overview see [l[ . Many schemes have been de- 
vised for which the error decreases at a much faster rate, 
compared to a straightforward approach of using each 
channel separately 0, HJ ■ This is analogous to the prob- 
lem of transmitting a reference frame [5], [H, [^, 0| • This 
increase in the rate of estimation is possible both with 
@, II II H Hi and without @,[|| the use of entanglement. 

As far as we are aware, no work has been done on the 
situation in which there are n non-identical channels. In 
this paper we investigate whether an increase in the rate 
of estimation is possible in this case. We look at the 
following two cases: (i) there are n unitary channels be- 
longing to SU(d) specified by the parameters 9 1 , . . . , 9 n , 
(ii) there are n not necessarily identical one-parameter 
unitary channels which depend on a common variable. 
In the former scheme an increase in the rate of estima- 
tion is not possible; in the latter scheme it is. 

To quantify the performances of estimation schemes 
we use the asymptotic limit of the mean square error 
E[{9 1 - 6 l )(03 - 9 j )]. Using the maximum likelihood 
estimator, as the number of measurements N — > oo, 
the mean square error is approximately (1/N)F^ 1 (9) @, 
where Fm(9) is the Fisher information matrix obtained 
from a single measurement M. The Fisher information 
for the parameter 9 = (9 , . . . , 9 P ) is the p x p matrix 
with entries 

F M {9) jk = J P (t 9) { —. j 

The Cramer-Rao inequality states that the mean square 
error of an unbiased estimator t(x) is less than or equal 
to the inverse of the Fisher information, 

m.s.e. e [t(x)]>F M (9)- 1 . (1) 
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It has been shown p3| that the SLD quantum informa- 
tion H(9) is an upper bound on the Fisher information, 
i.e. Fm{9) < H{6), The SLD quantum information for 
the parameter 9 = (0 1 , . . . , 9 P ) is defined as the matrix 
with entries 

where X J is any self-adjoint solution to the matrix equa- 
tion 

^H(/^+AV). (2) 

Given a state pg depending on some unknown parameter 
9, and a POVM {M m } we get a measurement outcome 
x. The quantum Cramer-Rao inequality states that the 
mean square error of an unbiased estimator t(x) of 9 is 
less than or equal to the inverse of the SLD quantum 
information, i.e. 

m.s.e. g [t(x)] > F -1 (0). (3) 

When p > 2, there exist families of states pg and ag for 
which neither H(pg) > H{ag) or H(pg) < H(ag) is true. 
To deal with this, we compare the traces ti{H(9)} of the 
SLD quantum information. 

Previous work has looked at how the error scales with 
the number of times U is used. We look at how the error 
scales with the number of input states used, as this is 
more convenient for us. 



I. 'INDEPENDENT' CHANNELS 

First we look at estimating n unitary channels from 
SU(d). The jth channel is specified by the parameter 
9 3 = . . . , 9 3 d 2_ l ). These channels are supposed 'in- 
dependent' in that there is no functional relationship be- 
tween 9 1 , 9 2 , . . . , 9 n ~ 1 and 9 n . Estimation of the channels 
is a parametric problem. 

We look at the mapping 
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We show in Appendix A that using a tensor product 
of maximally entangled input states po — &j = iP 3 mes is 
sufficient to maximize the trace of the SLD quantum 
information of the output states of This is not a 
necessary condition, as the 2n-partite entangled state 
1 j\fd • |ej) (£)■■•(£) |e,-) also attains the maximum SLD 
quantum information. However, since these states are 
significantly harder to produce, we are better off using 
maximally entangled states. The SLD quantum informa- 
tion using a tensor product of maximally entangled states 
is attainable [Tl| , so asymptotically the mean square er- 
ror is (l/AOff- 1 (0). 

For the mapping ((4]) , an optimal estimation procedure 
for n 'independent' unitary channels, in terms of tr{_ff}, 
is to estimate each one individually using a maximally 
entangled input state. The optimality of this procedure, 
in terms of tr{H(9)} has been shown by Ballester [111 ]. 

II. 'DEPENDENT' CHANNELS 

We look at the case where we have n not necessarily 
identical channels which depend on a common parameter 
9. These channels are of the form 

U S = ( o e ih(6) \ ---> U e = ( o e i/„(0) ) ' ( 5 ) 

where < 9 < t and Ti = C 2 . We can use each chan- 
nel N times. We impose the following conditions on the 
functions ff (a) fj(9) :Knt, (b) dfj{9)/d9 > 0, and 

(c) < J2 3 fA°) < ^ for a11 3 and °- We look at the 
mapping 

Po » ®7=i(C7g ® I R ) Po {U$ <8> I R ), (6) 

where po € S(® n (T-t ($> H R )); we denote by S(H) the set 
of states on H. We compare the SLD quantum infor- 
mation between (i) using a tensor product of maximally 
entangled states, (ii) using a 2n-partite entangled state. 
In this case using 2n-partite entanglement gives consid- 
erably larger SLD quantum information. 

Using the input state po = p x mes ® • • • <8> Pmes > where 
Pies = WuWul m = 1/V2(|00> + |11)), gives an 
SLD quantum information of Y^j=i {dfj/d9) 2 , which is 

attainable using a tensor product of the POVM {Mq = 

\i&)wl\M =i-iv>ix^|}- 

If we use the input state po = IV'oXV'oIj where |i/>o) = 
1/V2(|00..0) + |ll..l)), we get an SLD quantum informa- 
tion of (X™ df 3 /d9) 2 . As dfj/dB > for all j, this is 
considerably larger than the SLD quantum information 
using a tensor product of maximally entangled states. 

The phase (f> of the output state is in one-to-one cor- 
respondence with 9, because of conditions (b) and (c); 
in fact cj) = £" =1 /j(0). Usin g the POVM {M = 
|Vto><M-Wi = I-hfoKlfo|} we get p{0-<t>) = cos 2 (0/2). 
We perform © N times and each time use this POVM. 
From the measurement outcomes we get an estimate (f> 



of (f> by cos 2 (0/2) = no/N, where n is the number of 
times we get the outcome x = 0. Because of condition 
(c) on the functions fj(9), cos 2 (</>/2) is in one-to-one cor- 
respondence with <f> and hence 9. This POVM gives a 
Fisher information equal to the SLD quantum informa- 
tion. Hence as N — + oo we get a mean square error of 
l/(^(Ei=i dfi/dQf) compared to l/{NYT =M ] /d9f) 
using a tensor product of maximally entangled input 
states. 

It is known that when we have n identical simple uni- 
tary channels, we can obtain an increase in the rate of 
estimation without using entanglement [l| . A simple way 
is to use each of the n channels in sequence on a single in- 
put state. Consider the unitary channel Ug = Diag(l, e l9 ) 
, i.e. 

po ^ Ugpo U'f , Ul l = U e Ue---U e . (7) 

We repeat (|7J| N times using the input state p x — 
|^)(^l where \i(j x ) = 1/V2(|0) + |1», and the POVM 
{\ip x ) (ipxl,^- — \4'x){4'x\}- We call this the sequential 
scheme. As N — > oo, we obtain a mean square error that 
scales as l/(Nn 2 ) . We can use this sequential scheme 
with n non-identical channels, i.e. 

p Q ~{Ul>-U$)p Q (Ui~-U$)\ (8) 

where po € S(H). Using the input state p x — 
\ip x )(ip x \ where \ip x ) = 1/V2(|0) + |1», and the POVM 
IV'xX^I} we get p(O;0) = cos 2 (0/2). Per- 
forming ([8|) N times we get an estimate 4> and hence 0. 
We get the same Fisher information as ([6]) using a 2n- 
partite entangled state, and hence the same mean square 
error. 

Since multi-partite entanglement is difficult to create, 
we are better off using the sequential scheme ([5J . 

Conditions (b) and (c) for the functions fj (9) are very 
strict. We can still get a considerable increase in the 
rate of estimation without these conditions. If condition 
(c) does not hold, the phase 4> of the output state is not 
in one-to-one correspondence with p(0;<fi) = cos 2 (</>/2). 
If we modify condition (c) to < JV fj(9) < 2tt, we 
can still find </> but we shall need to perform an ex- 
tra measurement {\ip y )(ip y \^I — \i/>y)(ip y \}, where \ip y ) = 
1/V2(|0) + «|1)), a small number of times to determine 
the sign of cos(0/2). Then we can estimate 4> and hence 
9. If < J^j fj(@) — does not hold, the phase <f> of the 
output state is not in one-to-one correspondence with 9. 
We can get around this using a method similar to that of 
Zhengfeng et al [H, based on Rudolph and Grover Q- 

Zhengfeng et al [l[ looked at the case where fj(6) = 
2tt9 for all j. Their scheme involves first using a sin- 
gle channel n times to get an interval in which 9 almost 
certainly lives. Then they use two or three channels in 
sequence to 'amplify' 9. This is repeated n times until 
they get a narrower interval for 9. This is continued un- 
til all channels are being used simultaneously. The mean 
square error scales as (log N' /N') 2 , where N' is the total 
number of times Ug is used. 
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In our case, the situation is more complex, as the func- 
tions fj{6) are more general and not necessarily identical. 
The finer details of how we go about this and how the 
mean square error would behave, depend on the func- 
tions fj(0). We give a very brief overview of a possi- 
ble procedure. We start by using a single unitary Ux n 
times to get an interval for f%(6) and hence 9. Then we 
use UiU\ to get an estimate of f\{9) + /2(d) and hence 
a more accurate estimate of 9. We continue this pro- 
cess till we are using all the channels simultaneously and 
we have a very narrow interval for 9. We expect that 
asymptotically, the mean square error is approximately 
V(^(Epi<!f3/^) 2 )i where N is the number of input 
states used. We leave a more in-depth analysis for fur- 
ther work. 

If neither (b) nor (c) are satisfied we propose the fol- 
lowing scheme: (i) Use each of the channels individu- 
ally n times to get an estimate 9, (ii) Divide the chan- 
nels into two groups: A — {U J g .dfj/d9 > at 9}, 
B = {Ug,dfk/d9 < at 9}, (iii) Use an iterative pro- 
cedure for the two groups separately, but sharing infor- 
mation about 9 to make the confidence intervals shorter. 

APPENDIX A: PROOF 

We are looking at unitary channels of the form 
Ug = exp(iJ2j Qjtj) where tj = Vp tr{t,} = and 
trjtjifc} = Sjk- We denote by H S j (po), the SLD quantum 
information for the jth channel 

Po ' * (Uoj <8>I)A>(0flj 

The SLD quantum information of using a tensor prod- 
uct of maximally entangled states is Hg{®™ =1 p mes ) = 
Diag(7J e i(p mes ), . . .,Hen(p mes )). 

Lemma 1 For all unitary channels of the form 

Po (Ug <8> T)pa(Ug (g> 1), (Al) 

the trace of the SLD quantum information is maximized 
by a maximally entangled state, i.e. 

tr{Hg(p Q )} < tr{H g (p mes )}. (A2) 

Equality holds in IAS\) if and only if po is a maximally 
entangled state. 

Proof. From the appendix of Ballester [jjj we know that 
for the channel (|A1[) 

^{{Heipmes^Heipo)} <d 2 -l. (A3) 

It is simple to show that for unitary channels of the form 
Ug = ex-V(iJ2j@jtj) we have Hg(p mes ) = (4/d)I<p_i. 
Substituting into (|A3[) we get 

tr{Hg(p Q )} < [ d ' = tr{Hg(p mes )}. (A4) 



Since equality holds in (IA3|) if and only if po is a maxi- 
mally entangled state [12], equality holds in (|A4p if and 
only if po is a maximally entangled state. 

Lemma 2 The trace of the SLD quantum information 
for 01) is maximized by a tensor product of maximally 
entangled states , i.e. 

t?{Hg(po)} < tr{H e (®] =1 pl nes )}, (A5) 

Po G S(® n (H ® H R )), p 3 mes G S(H ® H R ). 

For pure states a solution of (|2|) is A J = 2dpgjd9K It is 
not difficult to show that for the set of states UgpoUg the 
SLD quantum information is the matrix with entries 

Hg jk (p Q ) = 4mr{[/> (7 e fct } 

+ Atr{U 3 e poUl}tr{UlpoUl}, (A6) 

where U 3 e = dg 3 Ug. In Q Ug = <E>] =1 U 9] <g> l R . 

Consider an arbitrary diagonal element of the SLD 
quantum information, i.e. Hg mm , where m(j,k) = (j — 
l)p + k, corresponding to the parameter 9\. From (|A6|) . 

He mm { P o) = 4tr(C/^po^ fct ) + 4 (tr(C/f p Q Ul)f , 

where W k = d g j Ug. For convenience we label the Hilbert 
space on which Ugi <8 Ir acts as Hi and so the Hilbert 
space on which Ug acts is <S>™ =1 Wi. Taking the partial 
trace over Kb, = we get 

Hg mm (Po) = 4tr((C^>(g)I)p A (C7« t (g)l)) 

+ 4(tr(^®I)p A (c7j,®l)) 2 , (A7) 

where p A = tr^ B .{po}, the reduced state of po- In gen- 
eral, p A is mixed, and so p A = Y^i=iPiPi: where pi are 

pure and orthogonal. Substituting p A = Yli=iPiPi m t° 
(E3),weget 

d2 ( 

i=l ^ 

+ 4(tr([/ Jfe ®I)pi{tf ®I)) 2 \ 
and hence by (IA6|1 

d 2 

H mm (Po) = ^PiiHgj {pi))kk- 
i—1 

Summing over k we get 

d 2 -l d 2 

k k=l i=l 

d 2 

= ^2pitr{H 0j (pi)}. 
i=l 

< tr{H B3 (p mes )}. 
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Summing over j we get (|A5[) . 
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